Test of a simple and flexible Ss model molecule in a-Sg crystals 
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a-Ss is the most stable crystalline form, at ambient pressure and temper- 
ature (STP), of elemental sulfur. In this paper we analyze the zero pressure 

Q low temperature part of the phase diagram of this crystal, in order to test 

xn \ 

a simple and flexible model molecule. The calculations consist in a series of 

^ ■ 

molecular dynamics (MD) simulations, performed in the constant pressure- 
£h ' constant temperature ensemble. Our calculations show that this model, that 

o 

gives good results for three crystalline phases at STP and T>300K, fails at low 
temperatures, predicting a structural phase transition at 200K where there 



o 

(N 



should be none. 



Introduction: 

Elemental sulfur is of relevance in geophysics, astrophysics, material sciences and mas- 
sively used in industry fflJ2]. Recently, a high pressure superconductor transition, with 



T C =17K, has been found ||. 



Elemental sulfur is widely known by the unusual characteristics of its liquid phase. Nu- 



^ ■ merous works have been dedicated to the study of its liquid-liquid (A) phase transition, due 
to a large variation of its molecular composition. The A transition is of relevance in numer- 



ous fields, volcanology is an example [|]-|6[ . The liquid phase is also of enormous theoretical 
interest, with studies ranging from models to describe the A transition 0^0? to molecular 
dynamics (MD) simulations of the liquid [|T^-[T^| and to calculations of the stability and 
structure of isolated S n molecules (T3|-[R|- 

In spite of the numerous experimental and theoretical studies on the liquid phase, the 
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phase diagram and properties of the crystalline phases of elemental sulfur are far from well 
known 0,0] ■ A simple inter- and intramolecular potential, able to reproduce the structural 
and dynamical properties of the crystalline phases, should be very useful to understand the 
complex experimental data of sulfur crystals. 

The most stable sulfur allotrope at ambient temperature and pressure (STP) is S 8 , a 
cyclic crown-shaped molecule, stable in solid, liquid and gas phases. Three crystalline phases 
have been clearly identified at zero pressure and T^400K, being orthorhombic a-Ss the 
most stable crystalline form at STP. Fig. la shows its "crankshaft" structure with four 
molecules in the primitive unit cell. The space group is (Fddd), with 16 molecules in 
this non-primitive orthorhombic cell. If a-3g is slowly heated, it shows a phase transition to 
monoclinic /3-S 8 at 369K, which melts at 393K. Nevertheless, most a-S 8 crystals do not 
convert to /3-Sg [fTTf] , they melt, instead, at 385. 8K fT]]^. 



Experimentally, the samples of pure a-Sg can be easily identified because these crystals 
are greenish- yellow at 293K and turn white below 193K. Commercial sulfur, instead, contains 
small amounts of S7 and conserves the yellow color at low temperatures J1J. The change 
of color is assumed to be due to a shift of the electronic band gap between valence and 
conduction bands . 

The firsts molecular dynamics (MD) simulations of a-Ss crystals were performed in the 
constant volume- constant energy ensemble, at normal pressure and at several temperatures 
fl8fl , with the only aim to study the Raman and ir. intramolecular frequencies in a crystalline 
enviroment. Numerous works ]nj have also been dedicated to the measurement and/or 



calculation, by lattice sums and in different approximations, of the k=0 vibrational modes. 

The first attempt to reproduce the zero pressure phase diagram of a-, [3- and 7-Ss crystals 
has been recently performed [pH|j , using a simple and flexible model molecule, based in that 



proposed by Venutti et al. p8[ . The model molecule reproduces the STP structures the 
three crystals, the measured heat of sublimation of a-Sg, and the orientationally disordered 
high temperature plastic phase of j3-S 8 crystals. 

In this paper we study the zero pressure and low temperature phase diagram of a-Ss crys- 



tals, as given by the simple and flexible model molecule of ref . |Tj| . The calculations consist 
of a series of constant pressure - constant temperature classical MD simulations, at several 
temperatures and zero pressure. This MD algorithm is useful to locate structural phase 
transitions, or in the search of the most stable structure for a given set of thermodynamic 
parameters P and T. 

Unfortunately, the model molecule of ref. [2(J fails to reproduce the measured low tem- 



perature structure of a-Sg crystals: a structural phase transition is found at 200K, were 
there should be none. Here we give a detailed account of the calculations performed at low 
temperatures on a-Sg crystals, of a disordered sample, and of the reasons for the failure of 
the proposed model. 



The inter- and intramolecular potential model: 

The main problem found in this type of calculations is the inter- and intramolecular po- 
tential model of this flexible molecule. The model of our simulations is a slight modification 
of that proposed in the constant volume MD simulations of a-Sg crystals, that succesfully re- 
produced the available experimental data on the crystalline intramolecular frequencies |TJ| . 



The intermolecular potential model of ref. jT3[] was of the Lennard- Jones (LJ) atom-atom 
type. In our constant pressure- constant temperature MD calculations, these LJ parameters 
|I"8| were slightly changed to improve the fit of the calculated STP crystalline structure and 
configurational energy of a-Sg to the experimental data. The parameters of our LJ poten- 
tial model, for non-bonded S-S interactions are: e= 1.70kJ/mol and cr=3.39 A. The cut-off 
radius of atom-atom interactions is 12 A and correction terms to the configurational energy 
and pressure, due to this finite value, are taken into account in the usual way: by integrating 
the contribution of an uniform distribution of atoms. 

A similar procedure is followed for the intramolecular potentials of the bending and 
torsional angles, that are based on those of ref. fllS] . Our flexible model molecule takes into 
account all the low frequency intramolecular modes that mix with lattice modes and can 
therefore be relevant to the onset of structural changes as a function of T and P values. 



Since the streching modes (from 400 to 500cm" 1 ) are always well above in energies than 
the rest of the vibrational modes (< 250cm" 1 ) and of the lattice modes (between 10 and 
~100cm _1 ), the S-S bond length of our model molecule is held constant at 2.0601A. The 
bending and torsional intramolecular motions (from 150 to 250cm" 1 ), instead, are taken into 
account and their mix with lattice modes is allowed in a straigthforward way. 
The intramolecular potential for bending angles S-S-S is taken harmonic 

V(p) = ±C fi {p-Po) a , 

with a force constant of C ( g=25700Kks/rad 2 and /3q= 108deg. The intramolecular potential 
for torsion (r) angles is a double well: 

V{t) = A t + B t cos(r) + C T cos 2 (r) + D T cos 3 (r), 

with A r =57.19 kg, B T =738.41 k B , CV=2297.90 k B and £> r =557.25 k B . These parameters 
describe a double well with minima at r—'t 98.8deg., and a height barrier of about 12kJ/mol. 



Calculations: 

The structural and dynamical properties of a -S 8 , are studied in the (N,P,T) ensemble, 
via a series of classical constant pressure - constant temperature MD simulations. The stan- 
dard MD algorithm allows volume and shape fluctuations of the MD sample in order to 
balance the applied isotropic external pressure with the internal stresses the temper- 



ature control of the sample follows the approach of Nose ||22| , |23| . The equations of motion 
of these flexible molecules are integrated using the Verlet algorithm for the atomic displace- 
ments and the Shake algorithm for the constant bond lengths constraints on each molecule 



24] , |25| . The final MD algorithm is identical to that used in a study of black Newton films 



The starting point of our simulations is the experimental structure of a-Ss at STP |27 
The sample was first equilibrated with a constant volume algorithm on a trajectory of 30000 
time steps (of O.Olps.). Afterwards, the runs in the (N,P,T) ensemble were performed by 



decreasing and increasing the temperature (up to 450K) in steps of 25 or 50K. At each point 
of the phase diagram the sample is equilibrated for 20000 to 30000 time steps and measured 
in the following lOps. Near the phase transitions the equilibration times were increased 
several times. 

The sample of a-Sg crystals consisted of 3x3x2 orthorhombic cells (288 molecules). Some 
points of the phase diagram were recalculated with a sample of 4x4x2 cells (512 molecules). 



Results: 

a) The a -phase: 

a-S 8 crystallizes, at STP, in the orthorhombic space group F d (D|£) with 16 molecules 
in a face-centered cell, the lattice parameters at 300K are a= 10.4646A, b=12.8660A and 
c=24.4860A ||27|| . The primitive cell contains 4 molecules at C2 sites. The molecular rings lie 
parallel to the crystallographic c axis and are alternately oriented parallel to (110) and (110) 
crystalline planes. Fig. la includes part of our sample at 300K and shows the structure in the 
[110] direction. At 300K our calculated lattice parameters 10.34(10)A, b=13.20(3)A 

and c=24.20(5)A, the largest difference is found in the value of b axis, 2 % higher than its 
experimental value. The molecules location and orientation follow the experimental data 
|27fl . At 300K the calculated value of configurational energy is -103.5(2)kJ/mol, from which 
the heat of sublimation can be estimated as lOl.OkJ/mol and compares well with a measured 
value of 101.8(2)kJ/mol f2§. 

Figs. 2a and 2b show the calculated vibrational density of states of the isolated molecule 
and that of a-Sg at 300K. They reproduce the experimental pattern, although the intramolec- 
ular modes are calculated about 10% higher than the experimental ones. 



Figs. 3a and 3b show the calculated configurational energy and volume per molecule as 
a function of temperature, for both samples of a-Ss- Unfortunately, our orthorhombic a-Sg 
sample of 288 molecules turns out to be unstable for T<200K and distorts to a monoclinic 
cell. The final molecular array closely resembles that of a-Sg (Fig. lb). Fig3c shows the 
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calculated angles of the unit cell as a function of temperature. This result was later on 
confirmed with a larger sample of 512 molecules (4x4x2 orthorhombic cells). Below 200K 
this large sample do not distort to monoclinic, but nevertheless its configurational energy 
and volume per molecule are higher than those for the sample of 288 molecules (Figs. 3a 
and 3b). The fluctuations of energy and volume, in the sample of 512 molecules, are also 
larger at T<200 than those at 200K. 

The coincidence of the experimental change of color at 193K and our calculated tran- 
sition at 200K, made us review the experimental data on this phase. The usual reference 
quoted for the lattice parameters values, as a function of temperature, is ref. |[30|| . Their 
reported values were determined by X-ray powder diffractrometry, on different samples and 
by measuring only three lines at 8 temperatures between 100 and 300K. Our calculated 
powder diffraction patterns, not included here, are quite similar for the orthorhombic and 
monoclinic samples and no conclusive answer can be obtained from these data. Neverthe- 
less, there is a careful single crystal determination of the structure of a-Sg at 300 and 100K 
ST|, where the crystalline charge distribution is measured. Both measured structures be- 



long to the same crystalline space group and the molecular arrangement is identical. This 
study clearly disregard any possible structural phase transition and shows that the simple 
model molecule, in spite of the amount of experimental data that it is able to reproduce for 
T>200K, should still be improved p0 . 



b) Phase transitions and the liquid phase: 

The experimental evidence shows that the phase transitions of Sg crystals are obtained 
when defects are present, if they are not, metastable states can be maintained for days 
p|J2]fr^] . This evidence is valid for the a — (3 transition, and also for the solid-liquid phase 
transition 0. In the last case, the disorder is generated because the S 8 molecules start to 
dissociate when the temperature is increased. 

In our simulations, all calculated ordered samples of a-Sg crystals melt at T>450K. That 
the disorder is essential to promote the solid-liquid phase transition can be checked by sim- 
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ulating a disordered sample of Ss molecules. A cubic sample of 216 disordered Sg molecules 
was studied in the range 350-450K, following the same procedure as for the other samples. 
In this case the MD algorithm allowed changes in the volume of the MD box, but not in its 
shape. Fig. 4 shows the calculated discontinuity in the configurational energy and volume 
per molecule at 400K, near that measured for a-S$ (386K) and /3-Sg (393K) crystals. For 
T>400K the sample is liquid, as was determined by its diffusion constant. These simulations 
should be valid for temperatures above and near the transition, when the fraction of broken 
molecules is small, and show that the model molecule is useful in this range of temperatures. 



The calculated distances between firsts neighbor atoms compares well with those deter- 
mmined from the pair distribution function g2(r) obtained by neutron measurements ||29|| . 
Fig. 4a includes our calculated bonded and non-bonded atom-atom distances of crystalline 
a-Ss at 300K, Fig. 4b corresponds to our disordered 'liquid' sample at 425K. The 3 th and 5 th 
peaks at 3.9A and 5.4A, which are not found in the liquid phase, correspond to non-bonded 
distances. Their disappearance, observed in the liquid phase, cannot be unequivocally at- 
tributed to broken bonds [^J. 



Conclusions: 

In this paper we study the zero pressure low temperature range of the phase diagram 
of S§ crystals, as given by a simple and flexible model molecule. The molecule model, a 



simplified and slight modification of that proposed in ref. ||18|| , is extremely simple and gives 
good account of many experimental facts for T>200K, including the experimental data on 
the structure and dynamics of a-, (3- and 7-Ss crystals around STP and the orientationally 
disordered plastic phase of /3-S 8 crystals [f2(]]. It also reproduces the solid-liquid phase tran- 
sition, near the experimental temperature, but only for a disordered sample of S§ molecules 
||20|| . The last calculation is only of theoretical interest, and we do not intended, here, to 
simulate the real liquid phase of elemental sulfur. Nevertheless, the pair correlation function 
calculated in this theoretical 'liquid' shows that the 3 th and 5 th neighbours identified in solid 



phases correspond to non-bonded atoms, and that their disappearance at high temperatures 
can be expected even if the molecules do not dissociate in short segments. 



In spite of all the properties of S 8 crystals that the simple and flexible model molecule 
reproduces, below 200K this model is unable to calculate the experimental orthorhombic 
structure of a-S$- Although no structural phase transitions are measured for a-S$ crystals, 
our sample of 288 molecules is unstable and transforms to a monoclinic structure. 
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Figures: 

Fig. 1: Calculated stable structures, as given by this model molecule: a) orthorhombic 
a-Ss at 300K, b) monoclinic a'-Sg at 150K. 

Fig. 2: Vibrational density of states: a) Isolated molecule and b) a-S s crystal, both at 
300K. 

Fig. 3: Calculated (a) configurational energy and (b) volume per molecule of a-Ss (large 
circles: sample of 288 molecules and small circles: sample of 512 molecules) and a disordered 
cubic sample (triangles: increasing T, squares: decreasing T) as a function of temperature. 
For T>400K, the disordered cubic sample is liquid, (c) Unit cell angles of the a-S 8 sample 
of 288 molecules. The lines are a guide to the eyes. 

Fig. 4: Pair correlation functions: a) a-S 8 at 300K, b) our theoretical 'liquid' sample of 
S 8 molecules at 425K. Full lines: total atom-atom g2(r) functions, dotted lines: excluding 
intramolecular distances. 
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